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• Nature of physical problem: Higher order leading logarithmic QED 
radiative corrections to W pair production and decay at high energy e + e~ 
colliders, including finite width of the W's and anomalous trilinear gauge 
boson couplings. 

• Method of solution: Monte Carlo event generation. 

• Restrictions on the complexity of the problem: The matrix ele- 
ments for the hard subprocess contain only the doubly and singly res- 
onant contributions constituting the so-called "CC11" diagrams, in the 
Born approximation. The hadronization of final state quarks is handled 
by external programs. Interfaces to the most popular QCD Monte Carlos 
are provided. 

• Typical running time: The time needed per unweighted event depends 
on the beam energy, cuts and CPU. The test run (10000 unweighted par- 
ton level events at LEP2 energies) takes approximately 420 CPU seconds 
on an aging 486DX2/66 PC running Linux 2-e. On contemporary personal 
computers, the 42 ms/event will be reduced to 5 ms/event and work sta- 
tions will need less than 1 ms/event. The CPU time per event increases 
mildly with the beam energy. 
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1 Introduction 



The spectacular success of the high energy electron positron colliders LEP at 
CERN and SLC at SLAC has confirmed the predictions of the Standard Model 
(SM) for the interactions between the gauge bosons and the fermions even at the 
level of electroweak radiative corrections. The high precision of the experimental 
data allows to put limits on parameters of yet unobserved particles like the 
Higgs boson through the appearance of these particles in weak loop corrections. 
On the other hand, the non-Abelian structure of the gauge sector of the SM 
predicts couplings between the electroweak gauge bosons which have not been 
tested directly. The forthcoming upgrade of LEP |, §] to LEP2 §, g with 
a center of mass energy of up to y/s « 200 GcV and a future 0.5-1.5 TeV 
e + e~ linear collider || (NLC, for short) will in particular make the trilinear 
WWj and WWZ couplings directly observable through their contribution to 
the production of W pairs. 

The main task of LEP2 will be a precise determination of the mass and the 
width of the W and the production cross section in the threshold region Q. 
At NLC, a measurement of the total cross section, the W angular distributions 
and in particular of their longitudinal polarization component will give a handle 
on possible anomalous couplings among the electroweak bosons induced by new 
physics beyond the SM ^ , and in turn reveal an insight into the mechanism 
of electroweak symmetry breaking. 

However, possible new physics is already severely constrained by present 
LEP and SLC data, and the effects to be expected at NLC (even more so at 
LEP2) are small. In order to extract these small effects, one has to have a 
precise knowledge of the radiative corrections within the SM. 

The electroweak radiative corrections to theproduction of on-shell W's to 
one-loop order are by now well establishe d P , |l0| . The influence of the finite 
width of the W's has been investigated in |12|. Also, the higher order QED 
corrections have been calculated in the leading log approximation (LLA) and 
their importance has been emphasized in ref. |13| . An exhaustive overview of 
the standard model predictions has been provided in [B and Q . 

Unfortunately, the experimental reconstruction of the W's and the deter- 
mination of their polarization is complicated by the fact that they may decay 
either into leptons with an escaping neutrino, or into hadrons, where the jet 
energies may be poorly known due to undetected particles. In addition, the 
radiative corrections due to emission of photons produce a systematic shift of 
the effective center of mass energy towards smaller values. Such effects may 
best be studied with the help of a Monte Carlo event generator. 

Although quite a few semianalytical calculations of the corrections men- 
tioned above have been available for some time, version 1.0 of WOPPER |l[| had 
been the first publicly available and complete Monte Carlo event generator. In 
the meantime, many new programs have been released [^J. 

This paper describes version 1.5 of the Monte Carlo event generator WOPPER. 
It is based on the lowest order cross section for the process e + e~ — > (W^+VK - ) — > 
4f and focusses on QED radiative corrections in the LLA resummed to all orders 
in a and the effects of finite width of the W's. The four- momenta of the exclusive 
hard photons are generated explicitly and treated with full kinematics. The W 
decays into fermions are implemented at the parton level, including leading 
QCD corrections to the W branching fractions. 
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This write-up is organized as follows: In section g we outline the physics 
underlying the algorithms implemented in WOPPER. The actual implementation 
is described in section ^[ The parameters controlling the execution of WOPPER 
are discussed in detail in section ^ and the FORTRAN-77 interface is presented in 
section ^[ Distribution notes, a listing of all external symbols and the output of 
a test run can be found in the appendices. 



2 Theoretical Background and General Features 

2.1 QED corrections at very high energies 

In the structure function formalism jl6| [l?], [ll| the expression for the radiatively 
corrected cross section reads 



er(s) = J dx + dx- D(x + ,Q 2 )D(x-,Q 2 ) a(x + x-s) , 



(1) 



where a is the Born level cross section of the hard process, D(x, Q 2 ) are the 
structure functions for initial state radiation, and Q 2 is the factorization scale. 

The structure functions D sum the numerically most important leading log- 
arithmic contributions 



a 



- log — re 6% (at LEP2 and NLC energies) 



(2) 



to the electromagnetic radiative corrections to all orders. They satisfy the QED 



evolution equation |16 



Q 



dQ' 



:D(x,Q 2 ) = 



1 



PrJz) 



l + z 2 

1 - z 



with initial condition 



D{x,m 2 e ) = 5(1 -x) 



(3) 
(4) 
(5) 



The solution to eq. (^|) automatically includes the very important exponentiation 
of the soft photon contributions to the radiative corrections. 

An explicitly regularized version of (|j|) used in the Monte Carlo implemen- 
tation is given by 



(0) 



a 
2^ 



dzP ee (z) 



D(x,Q 2 



It is crucial to note that the characteristics of the generated event sample do 
not depend on e, if it is chosen well below the experimental threshold for the 
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detection of soft photons. Theoretically, the cross sections will remain positive 
for all values of e, but limited storage for ultra soft photons and limited floating 
point range impose a lower limit. 

The radiatively corrected cross section (|lj) is implemented in a Monte Carlo 
event generator by solving the integro-differential equation ([j]) by iteration and 
taking into account the energy loss in the hard cross section. It is very similar 
to algorithms for quark fragmentation in QCD [jl9). As a by-product of the 
branching algorithm, the four-momenta of the radiated photons are generated 
explicitly. 

The initial state branching algorithm has already been used in the Monte 
Carlo generator KRDNOS [^o) and an improved version in the generator UNIBAB 
where the implementation has been described in detail. 

The momenta of the electron and positron after initial state branching are 
then used as the input momenta for the subgenerator of the hard subprocess 
described below. 



2.2 Born Cross Section 

In the general case, there are many Feynman diagrams contributing to the pro- 
cess e + e~ — > 4 fermions at high energies, even if one requires that the quantum 
numbers of the final state fermions be consistent with W pair production Q . 
However, the contribution of the individual Feynman diagrams may be easily 
estimated by counting the number of resonant propagators, where an interme- 
diate vector boson may come close to its mass shell. Using this naive estimate, 
one finds that the contribution of these 'background diagrams' is suppressed by 
a factor Tw/Mw ~2.5% for each non-resonant boson propagator, and may be 
reduced further by appropriate cuts on the invariant masses of the reconstructed 
W's. In fact, a full calculation |22] finds that this estimate is numerically jus- 
tified for LEP2 energies, unless electrons are in the final state and no invariant 
mass cuts are applied. 

In the current version 1.5 of WOPPER, we have taken into account the doubly 
and singly resonant Feynman diagrams shown in figures |^ and ^, also known 
as the "CC11" diagrams Q. Note that the combination of the fermions 1 and 
2 has quantum numbers consistent with a W + , while the combination of the 
fermions 3 and 4 has quantum numbers consistent with a W~ , even for the 
singly resonant diagrams. We will refer to these pairs as "pseudo-VF's" below. 
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e + u e + u 




d 



Figure 3: The "CC11" subset of Feynman diagrams contributing to four fermion 
production e + e~ — » qq'q"q"' . For qq'lvi and l~vil + vi final states only 10 and 
9 diagrams contribute, respectively, because of the uncharged neutrinos in the 
intermediate states. 

In the calculation of the Born amplitudes, we use an energy dependent width 
for both the W and Z propagators, that is given by the contributions of the 
light fermions to the vector boson self energies: 

, k 2 

ImE ff (fc 2 )E# -IV(fc 2 )M —2- M W T W , Re Y, w {k 2 ) ss (7) 

and similarly for Z's. 

It can easily be seen that, after integrating the modulus squared of the 
Born amplitudes over the decay angles of the pseudo W's, one obtains the same 
resonance formula for the total cross section for production of the pseudo W 
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pairs as in the well known doubly renonant case 11 , 23 , 24 , for each combination 
of final state fermions /1/2/3/4: 

a(12,34) = J ds+ds- ^ (7 off ( S ; S+ , S _;12,34) 

(8) 

Here a Q ^(s; s+, s_) denotes the off-shell cross section for the production of two 
pseudo-W's, and rV(s±) is the effective decay width of a virtual W into light 
fermions according to eq. (0) . The resonance factors are Breit-Wigner functions 
with energy-dependent width: 

1 1 

(9) 



D(s±) (s±-M^ + s ± T 2 w (s ± )- 

Note that equation (Q) also holds for the total cross section, after summation of 
over the decay channels of the pseudo W's: 

For the actual implementation it is useful to apply the mappings (7 = 
T w /Mw) 

e ± = arctan^ j (11) 



to eq. (10) in order to get a smooth integrand suitable for a Monte Carlo rejection 



algorithm: 

= f dt + d£-a(s; (12) 
2.3 Triple gauge bosons couplings 

In order to model the triple gauge boson vertices WWj and WWZ including 
possible deviations from the Standard Model ("anomalous couplings"), we use 
an effective Lagrangian as given by Hagiwara et al. HQ 

C W wv/gwwv = igj (Wl„W» -WlW%)V v + iK V WlW v V» v 
m w 

1 %^l^\y vX (13) 



rr> 



w 



Here V* stands for either the photon or the Z field, is the W field, W^ v 
d„W v - d v W^, = d^Vu - dvV^ and % v = \e^ Xa V x ° . 

ijiere we use the conventions of Itzykson and Zuber with e 0123 = +1. 
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These seven operators exhaust all possible Lorentz structures when we ne- 
glect the scalar components of the vector bosons, 

= , fl^W" = . (14) 

which is allowed if we couple the bosons to (almost) massless external fermions. 
In accordance with Hagiwara et al. || we choose 

gww~t = -e , gwwz = -e cot 9 W ■ (15) 

In the Standard Model at tree level we then have: 

Si — 9i — K ~« = K z = 1 ! ah others couplings = (16) 

We demand P and C conservation separately for the WW^f vertex, thus we 
have to set 

91 = 55 = ^7 = ^7 = • ( 17 ) 

Electromagnetic gauge invariance requires that gj is related to the charge of the 
W, 

gj = i , (is) 

so that we have only two free parameters (k 7 ,A 7 ) for the WW^ vertex, while 
there are seven for the WWZ vertex. 

For convenience, the anomalous couplings can also be specified in another 
popular parameterization (Sz, £7, xz, U-y, Hz, zz). The translation formulae are 
presented in section |4.2| below. 



2.4 Coulomb singularity 

Another class of universal corrections, which is important near threshold, is the 
so-called Coulomb singularity. We implement the correction |p5| 

ol-k A 2 f \[3 M + Af-F \\ 

crcouiomb = crBorn— T 1 arctan — (19) 

2/3 V 7T V 2plm/3M J J 

to the off shell total cross section (|To|). Stricly speaking, this correction should 
only be applied to the subset of doubly resonant diagrams. However, the 
Coulomb correction is only large for the W's close to their mass shell where 
the singly resonant diagrams are suppressed, and thus we apply it as a universal 
correction factor. 



2.5 Decays of the pseudo Ws 

In the case of the doubly resonant diagrams, the relative probabilities of the 
final state fermions are simply given by the branching fractions of the W's. 
This remains true even after taking into account the inclusive QCD corrections 
to the W decays to quarks to first order in as, because of the factorization of 
the amplitude into W production and decay. Neglecting fermion masses, one 
has: 

T(W -> qq') = N c (l + T(W -> lv t ) (20) 



However, in the case of the inclusion of the singly and doubly resonant diagrams, 
it is not clear how to take the QCD corrections into account, since no full calcu- 
lation is as yet available. Neglecting QCD corrections to final states with quarks 
is clearly unacceptable, and simply multiplying the contribution of the doubly 
resonant diagrams with the QCD correction factor spoils the gauge cancella- 
tions between singly and doubly resonant diagrams. Therefore, we have taken 
the approach to implement the QCD corrections in (|§|) by using an "effective 
number of colors" 

Nf = N c (l + ^) , 

which we refer to as the "naive QCD corrections" . 

The final state fermions of the decay W — > //' are therefore chosen with 
a probability according to the cross sections obtained from (|^), with the naive 
QCD corrections applied as described above. The branching fractions for hadronic 
decays into u, d, c, s, and b quarks are given by the corresponding CKM matrix 
elements, while decays into t quarks are assumed to be kinematically forbidden 
because of \m t — Mw\ S> ^w- Currently, the fermion masses are only taken 
into account kinematically and not in the matrix elements. 

The angular distributions of the pseudo Ws and of the decay fermions is 
generated by a standard mapping and rejection algorithm from the differential 
cross section da(s; S+, s_; 8; 12, 34)/cK7 and from the modulus squared matrix 
element, respectively. 

2.6 Hadronization 

If one or both of the VPs have decayed into quarks, they can optionally be 
hadronized using either the LUND string model p6| , the HERWIG cluster 
model |2^J or the ARIADNE J2^] color dipole model. In the present version 
it is not possible^] to study color-rearrangement effects in purely hadronic de- 
cays. Both W decays are handled separately. 

The hadronization model can be switched at run time with the parameter 
qcdmc, which takes the values 0, 1, 2 and 3. These correspond to no parton 
level, JETSET m, HERWIG Q and ARIADNE Q respectively. 

3 Implementation of WOPPER 1.5 

Like almost all Monte Carlo event generators, WOPPER is divided into three 
parts: initialization, event generation, and termination. These are described in 
this section. For ease of use, WOPPER comes with two application interfaces, so 
that a direct call to the lower level parts in this section will never be necessary. 
These interfaces will be explained in section 0. 

3.1 Initialization 

The initializations in WOPPER are used for computing the value of variables that 
will be used frequently during event generation. Examples are the calculation of 
electroweak couplings, the maximum of the off-shell Born cross section and other 

2 Unless one makes the necessary changes to the WOPPER sources in routines wwlund ( ) , 
wwhwigO andwwariaC). 



9 



internal steering parameters from the input parameters. This is accomplished 
by a call to the subroutine wwinit after setting the Monte Carlo parameters. 
Since WOPPER does not yet include weak corrections, most of the initializations 
performed are quite trivial. 

Finally, a standard /hepevt/ initialization record [Q is written, which may 
be used by the analysis program. 



3.2 Event Generation 

The routine wwgen produces an event on every call. The four momenta of all 
generated particles as well as supplemental information is written to a standard 
/hepevt/ event record, where it can be read by user supplied analyzers. See 



section 3.4 below for details on WOPPER's extensive use of /hepevt/. 

The first step is the generation of the initial state radiation by the branching 
routine wwbini using the algorithm of section and of the virtual masses of the 
pseudo W's occurring in the intermediate state. According to the cross section 
(fl2|), the rejection weight is calculated from the ratio of the actual off-shell cross 
section to the maximum determined in the initialization step, and the event is 
accepted with a probability according to this weight. 

After the effective center of mass energy and the pseudo W masses have 
been fixed, the quantum numbers of the event (helicity of incoming fermions 
and type of final state particles) are determined with probabilities proportional 
to their relative cross sections. The angular distribution and four-momenta of 
the pseudo W y s are generated in wwgww and wwgppr, respectively. Finally, the 
decay of the intermediate W's into the final state fermions is accomplished in 
the subroutine wwgdec. 



3.3 Termination 

The cross section for the generated events is obtained in the subroutine wwclos 
from the standard formula 

, & of successful trials , . 

a tot (s) = m^_{a(s;^)} ■ — - - — (21) 

where a is the off-shell cross section from eq. (|l2"|). s' varies between and s, 
and £± in the range allowed by the specified cuts. The corresponding statistical 
error from the Monte Carlo integration is given by 

Aa tot (s) = max {£(«';£+,£_)}■ (22) 
»';£+,£- 

/ (total # of trials — # of successful trials) ■ # of successful trials 
y (total # of trials) 3 

The cross section and the error on the cross section are placed into /hepevt/, 
where they may be read by the user-supplied analyzer. 



3.4 Additional information in /hepevt/ 

Because WOPPER uses the standard /hepevt/ event record internally, not only 
stable particles with isthep(i) = 1 will be present. Adapting the conventions 
of the HERWIG Monte Carlo [^Tj, we use the following status codes 
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• 101: e~ beam (positive z-direction), 

• 102: e + beam (negative z-direction), 

• 103: center of mass system of the collider, 

• 110: e + e~ hard scattering center of mass system, 

• 111: e~ before hard scattering, 

• 112: e + before hard scattering, 

• 113: virtual W~ after hard scattering, 

• 114: virtual W + after hard scattering. 

However, these entries have no physical significance and should never be used 
in any analysis (an exception to this rule are the beam particles 101 and 102, 
which are convenient for defining the reference frame and are used e.g. by the 
analyzer hepawk pS|| for this purpose). Only the particles with status code 1 
belong to the final state as predicted by WOPPER. 

If no hadronization Monte Carlo is active, final state quarks will be entered 
as "stable" particles with status code 1. 



4 Parameters 

The parameters controlling WOPPER version 1.5 are summarized in tables 
and 0. They will be described in the following subsections. 



4.1 Electroweak Parameters 

Since the present version of WOPPER does not include an electroweak library, the 
electroweak parameters, namely the masses of the electroweak bosons (masslw, 
masslz), their widths (gammlw, gammlz) and the Weinberg angle (sin2w) are 
treated as independent parameters. They enter an effective Born cross section 
and may be set directly by the user. As a special case, setting gammlw to re- 
produces the narrow- width approximation with on-shell W 's in the intermediate 
state. 

The bulk of the non-electromagnetic radiative corrections can be incorpo- 
rated into the hard cross section by using the running QED coupling cxqed (4M^-) ~ 
1/128 at the W scale. This value is, however, not correct for the initial state 
radiation of on-shell photons, where C(qed(0) ~ 1/137 has to be used. The 
inverse of the former value can be changed with ahpla and that of the latter 
with ahplaO. 

The presence of initial state radiation can be toggled using the parameter 
bstyle. The supported supported values are and 1, corresponding to no QED 
radiative corrections and LLA resummed initial state QED radiative corrections. 

The variable cell controls the number of Feynman diagrams taken into 
account. The default value is .true., which amounts to using all diagrams 
shown in figures || and |[ Setting this variable to false restricts the calculation 
to only the doubly resonant ("signal") diagrams in figure also known as 
"CC03" diagrams. 
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Variable name 


semantics 


Default value 


ebeam 


e + /e~ beam energy 


250 GeV 


epol 


longitudinal e~ beam polarization 





ppol 


longitudinal e + beam polarization 





scheme 


renormalization scheme 


1 


ahpla 




128 


ahplaO 


1/o:qed{0) 


137.0359895 


massle 


m e ± 


0.51099906- 10" 3 GeV 


masslw 


M w ± 


80.26 GeV 


gammlw 


Tw 


2.08 GeV 


masslz 


M z o 


91.1884 GeV 


gammlz 


Tz 


2.492 GeV 


gf ermi 


G F 


1.16639- 10- 5 GeV 


sin2w 


sin 2 9w 


0.2310 


alphas 


a QCD (M^) 


0.12 


cell 


all CC11 diagrams 


. true . 


f loops 


fermion loop contributions 


.false . 


coulom 


Include final state coulomb corrections 


.false . 


ckmvus 


V us = sin 9c 


0.2196 


ckmvcb 


v cb 


0.0400 


c km vub 


v ub 


0.0032 



Table 1: Standard model parameters controlling WOPPER. 



As a matter of convention, all gauge boson widths are treated in WOPPER as 
running widths (c.f. eq. (0)). However, it is well known that the inclusion of an 
energy-dependence in the width of the gauge bosons distorts the cancellations 
between the photon and Z-exchange contributions of the signal diagrams. It has 
been noted p0[ that the cancellations may be restored by including a minimal set 
of higher-order diagrams, namely the fermion loop contributions to the ZWW 
vertex. Neglecting the masses of all fermions except for the top-quark which 
is assumed to be very heavy, this so-called "fermion loop scheme" amounts to 
replacing the Standard Model ZWW tree-level vertex by: 

r ZWW _^ r ZWW _ r ZWW ( -i , /no\ 
1 ^J^floops — -Ltrco ' \ Mz ) 

The inclusion of this factor is controlled by the variable f loops. 

Since version 1.3 it is possible to choose canonical input parameters for 
benchmarking LEP2 Monte Carlos Q by changing the value of scheme. The 
following values are supported: 

• scheme = 0: "free scheme", all parameters are taken from the input and 
treated as independent parameters. 
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Variable name 


semantics 


Default value 


acpara 


parameterization 





kappag 




1 


lamdag 


A-y 





kappaz 


«Z 


1 


lamdaz 


Az 





glz 


g? 


1 


g4z 


flf 





g5z 


flf 





kapzt 


kz 





lamzt 


Az 





acdz 


<5z 





acxg 







acxz 


xz 





acyg 







acyz 


yz 





aczz 


zz 






Table 2: Anomalous couplings parameters controlling WDPPER. 



scheme = 1 (default): "Gf scheme" 



sin 2 ^ = ™^ M w) (24) 
IV = ^k( 3+ 2 -^CR) (25) 



V8vr 



scheme = 2: "sin#^- scheme" 



G F = "-pEE^MM (26) 
V2sin 2 6» w M^ 

rv = ^Mk( 3+ 2 -^£R) (27) 

V87T V 77 J 

scheme = 3: "Born scheme", tree level formulas, independent sin 6w, Gf, 
a QED (4M^) = ckqed (0) and 



w 



GpM^ ( 2a QCD 



V8ir 



3 + ^ (28) 



scheme = -1, -2, -3: same as the positive values, except for FV which 
is not calculated from the standard model expression but taken from 
gammlw instead. 
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Variable name 


semantics 


Default value 


cutmin 


minimum W virtuality 


GeV 


cutmax 


maximum W virtuality 


-s/S = 2EBeam 


nevent 


Number of events 


10000 


cc 


apply canonical cuts 


.false. 


cclvl 


level of canonical cuts 





ccver 


version of canonical cuts 


1 


bstyle 


Key for QED radiation 


1 


circe 


beamstrahlung 


. false . 


ciracc 


accelerator 


2 (TESLA) 


cirver 


version 


1 


cirrev 


revision 


19960401 


circht 


chattines 


1 


epsiln 


Internal infrared cutoff 


10~ 5 


qcdmc 


Key for QCD Monte Carlo 





rangen 


Random number generator 


1 


rseed 


Random number seed 


54217137 


rlux 


'Luxury' of generator 


3 


errmax 


maximum error count 


100 


verbos 


verbosity 





run id 


run identification 




stdin 


standard input 


5 


stdout 


standard output 


6 


stderr 


standard error 


6 



Table 3: Technical parameters controlling WDPPER. 



4.2 Anomalous couplings 

The parameterization of anomalous WWZ and WWZ couplings in table || is 
controlled by the acpara variable: 

• acpara = 0: no anomalous couplings at all. During initialization, k 7 , A 7 , 
K z, ^z, 9i , gf ! gi i Hz, will be reset to their standard model values. 

• acpara = 1: anomalous couplings in the Hagiwara et al. Q parameteri- 
zation, using the values of k 7 , A 7 , kz, Xz, 9i , gf ', gi , Xz specified in 
the variables from table |[ 

• acpara = 2: anomalous couplings in the (5z,x jl xz,y'y,yz,zz) parame- 
terization. The corresponding Hagiwara et al. parameters 

k 7 = 1 + x 7 (29) 
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A 7 = y 7 (30) 

n z = 1 + (x z + S z ) tan9 w (31) 

^z = y z tan9w (32) 

grf = l + <5 z tan(V (33) 

flf = (34) 

g§ = -— r z z ta,n6 w (35) 

H z = (36) 

A z = (37) 

will be used in the Monte Carlo. 



The relation (f35|) between and zz is only approximate. z z is defined as the 
coefficient in the higher dimensional operator 

Cpy = ^-d^p {d v w~ p w +p ~ d v w-»w +p + d v w +p w~ p - d^w+^w-p) 

M w 

(38) 

In the approximation of massless final state fermions and equal invariant mass 
of the intermediate Ws, the relation (^) can be shown. Radiative corrections 
cause a further deviation, because of the change in s. Until future version of 
WOPPER will support zz directly, ( |35|) is an excellent approximation. 

It should be noted that even a small value of z z will have a strong effect at 
NLC energies, because of the s-factor in ([35]) induced by the additional mass 
dimensions in Hgf- The natural value of z z will therefore be much smaller than 
that of the other parameters. 



4.3 Cuts 

In the present version 1.5 of WOPPER, only cuts in the virtualities of the inter- 
mediate Ws are implemented in the event generation: 

• cutmin: minimum virtuality of the intermediate W^s, 

• cutmax: maximum virtuality of the intermediate W^s. 
The cuts in virtualities have to satisfy the following conditions: 

< cutmin < cutmax < \/s = 2EBeam (39) 

A value of for cutmax will automatically be reset to the available maximum, 
namely 2E Be am- 

Since version 1.3 it is possible to apply the canonical cuts for benchmarking 
LEP2 Monte Carlos [Q to the event record by setting cc to . true . . These cuts 
will be reflected in the caclulated total cross section. 



4.4 Beamstrahlung 

WOPPER can be linked with the CIRCE library p] to include the effects of the 
beam-beam interaction at a linear collider. If the parameter circe is set to 
. true . , the momenta of the incoming electrons and positrons will be distributed 
according to the version (cirver) and revision (cirver) of the CIRCE parame- 
terization of a particular accelerator design class(ciracc). 
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4.5 Monte Carlo Parameters 



The remaining, more technical Monte Carlo parameters should be almost self 
explaining. Since our branching algorithm automatically includes soft photon 
exponentiation, the results will not depend on the value of the internal infrared 
cutoff epsiln (which is measured in units of the beam energy), provided it is 
kept well below the experimental energy resolution. However, it is not advisable 
to set it many orders of magnitude lower than the default value, because this 
may result in too high photon multiplicities that will overflow internal tables. 

A note on the random number generators available: the default value 1 of 
the parameter rangen corresponds to the standard RANMAR generator, that 
has been the generator of choice for quite some time. Recently, the quality of 
the random numbers generated by RANMAR has been questioned and unwanted 
correlations have been found, that caused large systematic errors in solid state 
physics simulations (3^]. A superior variation RANLUX has been proposed pifl , 
which is however much slower. If WOPPER has been linked with the CERN library, 
setting rangen to 2 switches to RANLUX which can be used at the 5 "luxury levels" 
to 4. At the highest "luxury levels" , RANLUX will be rather slow, however. 

However, since event generation involves a lot of decisions that effectively 
randomize the subsequences used by themselves, we do not expect that the cor- 
relations in RANMAR have any significant effect on the event samples generated 
by WOPPER. The RANLUX option has been added to WOPPER for some experimen- 
tation only. It is left in only because there is no particular reason for throwing 
it out again. 

An even more interesting option is the random number generator recently 
proposed by Donald Knuth |55|. This generator is selected by setting rangen 
to 3. It generates 30-bit integers with the following desirable properties: 

• the bit stream passed all the tests from George Marsaglia's "diehard" suite 
of tests for random number generators |36| (the result of the "birthday- 
spacing" test has to used with great care, because the position of the least 
significant bit of the 30-bit words varies if the bit stream in distributed 
into 32-bit words). 

• it can be implemented with portable signed 32-bit arithmetic (Fortran 
can't do unsigned arithmetic). 

• it can create at least 2 30 — 2 independent sequences, which are selected by 
a 30-bit integer seed value. 

• it is very fast. 

The default value of rangen is still 1, but all users are encouraged to give the 
Knuth generator a try. 

4.6 QCD Parameters 

The parameters for JETSET should be accessed trough WOPPER's lugive com- 
mand, which is translated directly to JETSET's LUGIVE () subroutine. See the 
JETSET manual p6| for a comprehensive description of the available parameters 
and their effects. 
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WOPPER 


HERWIG 


semantics 


Default value 


hwqcdl 


QCDLAM 


AqcD/GeV 


0.18 


hwrms 1 


RMASS(l) 


m d / GeV 


0.32 


hwrms 2 


RMASS (2) 


m u / GeV 


0.32 


hwrms 3 


RMASS (3) 


m s /GcV 


0.5 


liwrms4 


RMASS (4) 


m c / GeV 


1.8 


hwrms 5 


RMASS (5) 


m 6 /GeV 


5.2 


liwrms6 


RMASS (6) 


mtj GeV 


100.00 


hwrms 


RMASS (13) 


m e / f -/GeV 


0.75 


hwvqcu 


VQCUT 


Quark virtuality cutoff 


0.48 


hwvgcu 


VGCUT 


Gluon virtuality cutoff 


0.10 


hwvpcu 


VPCUT 


Photon virtuality cutoff 


-1.00 


hwclma 


CLMAX 


Max. cluster mass parameter 


3.35 


hwpspl 


PSPLT 


Split cluster parameter 


1.00 


hwqdiq 


QDIQK 


Max. scale for g — > diquark 


0.00 


hwpdiq 


PDIQK 


g — > diquark rate parameter 


5.00 


hwqspa 


QSPAC 


Spacelike evolution cutoff 


2.50 


hwptrm 


PTRMS 


Intrinsic px 


0.00 



Table 4: HERWIG parameters accessible from WOPPER. 



Since HERWIG does not sport the equivalent of the LUGIVEO routine, its pa- 
rameters have to be accessed through the standard WOPPER access mechanisms. 
Tables |] and || provide a list of the available HERWIG parameters and the names 
under whith they are known to WOPPER. See the HERWIG manual |27j for a com- 
prehensive description of the effects of these parameters. 

ARIADNE doesn't supply a high level parameter changing interface (like LUGIVE ( ) ) 
as well. Since the tunable parameters are more copious than in the case of 
HERWIG, we have refrained from making them accessible through the WOPPER 
parameter interface. Users who wish to change ARIADNE paramters are invited 
to add a call to their parameter changing routine to wwariaO. The standard 
ARIADNE tuning sets can however be accessed with the artune command, which 
passes its argument to the ARTUNE () routine. 

5 FORTRAN-77 Interface 

WOPPER version 1.5 provides two application program interfaces on different lev- 
els. The higher (preferred) level consists of the command interpreter wwdcmd 
that accepts commands in form of character* (*) strings. This driver commu- 
nicates with the analyzer hepawk (2^] by default. The lower level consists of two 
FORTRAN-77 subroutine calls: wwpsrv and wopper. 
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WOPPER 


HERWIG 


semantics 


Default value 


hwipri 


IPRINT 


Printout option 


1 


hwmaxp 


MAXPR 


Max. printouts 





hwmaxe 


MAXER 


Max. errors 


10 


hwlwev 


LWEVT 


Event output unit 





hwlrsu 


LRSUD 


Sudakov input unit 





hwlwsu 


LWSUD 


Sudakov output unit 


77 


hwsudo 


SUDORD 


Sudakov order in as 


1 


hwnrnl 


NRN(l) 


First random seed 


17673 


hwnrn2 


NRN(2) 


Second random seed 


63565 


hwazso 


AZSOFT 


Soft gluon azimuthal correlations 


.true. 


hwazsp 


AZSPIN 


Gluon spin azimuthal correlations 


.true. 


hwblli 


B1LIM 


i?-clustcr — > 1 hadron parameter 


0.0 



Table 5: Technical HERWIG parameters accessible from WOPPER. 



* wopperappl.f 



call wwdcmd ('set ebeam 95') 
call wwdcmd ( ' set acpara 2 ' ) 
call wwdcmd ('set acdz 0.01') 



select y/s= 190 GeV 
select anomalous coupling 
5 Z = 0.01 



call wwdcmd ( ' init ' ) 

call wwdcmd ('generate 10000') 



initialize the generator 
generate 10000 events 



call wwdcmd ('close') 



cleanup 



Figure 4: Higher level FORTRAN-77 interface 



5.1 Higher Level Interface 

The simple commands understood by wwdcmd are (here keywords are typeset in 
typewriter font and variables in italics; vertical bars denote alternatives) 

• initialize 

Force initialization of WOPPER and write an initialization record into the 
/hepevt/ event record, which should trigger the necessary initializations 
in the analyzer. 

• generate [n] 

Generate nevent events and call hepawk to analyze them. If the optional 
parameter n is supplied, nevent is set to its value. 

• close 
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Write a termination record to /hepevt/, which should trigger the neces- 
sary cleanups in the analyzer. 



• statistics 

Print performance statistics (this is usually only useful for the WDPPER 
developers, who are tuning internal parameters). 

• quit 

Terminate WDPPER without writing a termination record. 

• exit|bye 

Write a termination record and terminate WOPPER. 

• set variable ival\rval 

Set physical or internal parameters. See the tables |l| and [}] for a compre- 
hensive listing of all variables. For example, the command ' set ahpla 127.0' 
will set the QED fine structure constant to 1/127. 

• print variable\&ll 

Print the value of physical or internal variables, print * causes a listing 
of all variables known to WOPPER and print prefix* causes a listing of 
all variables starting with the given prefix. 

• debug|nodebug flag 
Toggle debugging flags. 

• testran 

Test the portability of the random number generator. We use a generator 
of the Marsaglia-Zaman variety |3^ |, which should give identical results 
on almost all machines. 

• banner 

Print a string identifying this version of WOPPER. 

• echo message 

Print message on standard output. 

• lugive string 

Pass string to JETSET's LUGIVE () routine. See Q for details. 

• artune string 

Pass string to ARIADNE's ARTUNE () routine. See Q for details. 

Unique abbreviations of the keywords are accepted, i.e. ' g 1000' generates 
1000 events. The tokens are separated by blanks. Blank lines and lines starting 
with a # are ignored and may be used for comments. For portability, only the 
first 72 characters of each line are considered. 

On UNIX systems WOPPER reads the default startup files .wopper in the 
user's home directory and the current directory, if they exist. 
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* 



wwdcmd . f 

subroutine wwdcmd (cmdlin) 
character* (*) cmdlin 



else if (cmdlin. eq. 'init') 
call wopper (0) 
call hepawk ('scan') 



initialize WOPPER 
print initialization record 



else if (cmdlin. eq. 'generate') 
do 10 n = 1, nevent 



call wopper (1) 
call hepawk ('scan') 



generate an event 
analyze the event 



10 



continue 



else if (cmdlin. eq. 'close') 
call wopper (2) 
call hepawk ('scan') 



get total cross section, etc. 

finalize analysis 



else 



end 



Figure 5: Event generation loop 



5.2 Lower Level Interface 

The subroutine wopper (code) has a single integer parameter. The parameter 
code is interpreted as follows: 

• 0: initialize the generator and write an initialization record to /hepevt/. 

• 1: generate an event and store it in /hepevt/. If WOPPER has not been 
initialized yet, the necessary initializations are performed, but no initial- 
ization record is written. 

• 2: perform final calculations and write the results to /hepevt/. 

Figure ^ displays excerpts from a simplified version of wwdcmd that make the 
correspondences between the two levels of the FORTRAN-77 interface explicit. 

WOPPER's parameters can be accessed on the lower level by the subroutine 
wwpsrv (result, action, name, type, ival, rval, dval, lval). The pa- 
rameter is specified by its (lowercase) name in the character* (*) string name. 
The string action is either 'read' or 'write' corresponding to whether the 
parameter is to be inspected or modified. The type of the parameter ( ' int ' , 
'real', 'dble', or 'lgcl') is specified in type; it is an input parameter for 
'read' and an output parameter for 'write ' . Depending on this type the value 
is passed in ival, rval, dval, or lval, respectively. The following error codes 
will be returned in the string result: ' ': no error, 'enoarg': invalid action, 
'enoent': no such parameter, 'enoperm': permission denied, and 'enotype': 
invalid type. 

The protection scheme implemented with this parameter handling has been 
described in po|| . Its main purpose is to guarantee consistency of user defined 
and computed parameters in the generation phase of the Monte Carlo. 
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usage: $0 [options] [var=val] 



var=val 



assign VALue to WQPPER VARiable 



beam-energy energy: beam ENERGY (= sqrt(s)/2) in GeV 
sqrt(s) (= 2 beam energy) in GeV 
generate NUM events 

write a summary record every NUM events 
use hepawk SCRIPT 



-e energy 

-s roots, — roots roots: 
-n num, — events num: 
-r num, — report num: 
-H script , — hepawk script : 



-q yes I no, — QED yes I no: 
-Q mc, — QCD mc: 
-a mode, — ac mode: 
-f file, —file file: 
-R rng, — rangen rng: 
-S seed, — seed seed: 



apply QED raditiative corrections 
use QCD Monte Carlo MC 
select anomalous couplings MODE 
source FILE with variable settings 
use Random Number Generator 
random number SEED 



-k, — keep: 

-o file, — output file: 
-p file, — preserve file: 
-x cmd, — sed cmd: 
-t dir, — tmpdir dir: 
-b bin, — binary bin: 



keep the temporary run directory 
write output to FILE (implies -p FILE) 
salvage FILE from temporary run directory 
filter hepawk script through sed CoManD 
use DIR as temporary run directory 
use BIN as WOPPER binary 



-h, — help: 



this message 



Figure 6: Command line options of the wrun script for WOPPER on UNIX systems. 

5.3 Command Line Interface for UNIX Systems 

The distributed version of the WOPPER main program reads commands from the 
standard input channel and passes them directly to wwdcmd. WOPPER writes 
results to the standard output channel, user supplied analysis routines may 
write additional files with histograms, etc. 

Using no other operating system services than the standard input channel 
for input makes the program very portable. At the same time, it is not the 
most user friendly approach for repetitive tasks, like scanning parameter sets. 
Therefore we provide a shell script wrun, which can be used to control wopper 
from the UNIX commandline. 

In figure |^ the usage message of wrun is shown. Variable can be set on the 
command line. The wrun script will create an approriate input file for WOPPER: 

(var)= (value): sets the WOPPER variable (var) to the value (value). Ex- 
amples: ccll=false, sin2w=0 . 23103, etc. 

Some of these variables can be changed by an option as well: 

-e (energy) or — beam-energy (energy): set the beam energy, this is 
equivalent to ebeam=i?B- 

-n (num) or — events (num): set the number of events, this is equivalent 
to nevent=n. 
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-r (num) or — report (num): write a summary record every everz (ram) 
events, this is equivalent to report=rt. 

-S (seed) or — seed (seed): set the seed of the random number generator, 
this is equivalent to rseed=n. 

For some variables symbolic names or convenience translations are available: 

-s (roots) or — roots (roots): set the center of mass energy -y/s, this is 
equivalent to ebeam=Y / s/2. 

-q yes|no or — QED yes|no: toggle electromagnetic radiative correc- 
tions, — q yes is equivalent to bstyle=l and -q no is equivalent to 
bstyle=0. 

-Q (me) or — QCD (me): select a hadronization Monte Carlo: (me) can 
be parton, jetset, herwig or ariadne. 

-a (mode) or — ac (mode): select an anomalous couplings parameteriza- 
tion: (mode) can be no or SM for no anomalous couplings, Hagiwara, 
kappa, lambda for the k-X parameterization or delta for the S-x-y 
parameterization. 

-R (rng) or — rangen (rag): select a random number generator: (rag) 
can be ranmar, ranlux or Knuth. 

The behavious of the script can be changed with the following options: 

-b (bin) or — binary (bin): select a different WDPPER binary. 

-o (file) or — output (file): write the WDPPER output to (file) 

-t (dir) or — tmpdir (dir): specify a temporary run directory. By de- 
fault, WOPPER will be run in a temporary directory with a unique name. 

-p (file) or — preserve (file): copy the file (file) to the starting directory 
before deleting the temporary run directory. 

-k or — keep: don't delete the temporary run directory. 

The following options are useful if hepawk p9f is used for event analysis: 

-H (script) or — hepawk (script): select a hepawk script. 

-x (cmd) or — sed (cmd): filter the hepawk script through the UNIX 
stream editor sed. This is particularly useful for changing the name of 
the file in which hepawk saves its histograms. For example, if the script 
contains the command save ( "@@@") ; , the option -x s/@@@/f oo .hbook/ 
will cause hepawk to save the histograms in foo. hbook. 

This shell script can be used to run test suites easily. For example, the 
tests in can be run with the shell script in figure [?] provided that the file 
ac. hepawk contains an appropriate hepawk script. 
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#!/bin/sh 
run () { 

. /wrun -a delta ac$l -q $2 -e 95.0 -n 100000 ccll=false \ 
scheme=-2 ahpla=128.07 sin2w=0 . 23103 mass lz=9 1 . 1888 \ 
masslw=80.23 gammlz=2 . 4974 gammlw=2 . 08 \ 
-o ac.$3.out -p ac.$3.hbook -H ac.hepawk \ 
-x s/ac .hbook/ac . $3 .hbook/ 

> 

for q in 1; do 

run acpara=0 $q sm . q$q 

for d in 2 5; do 

for v in dz xg xz yg yz zz; do 
for s in - +; do 

run $v=${s}0.$d $q $v$s$d.q$q 
done 
done 
done 
done 
exit 

Figure 7: Running tests for various anomalous couplings with WOPPER on UNIX 
systems. 

6 Conclusions 

We have presented version 1.5 of the Monte Carlo event generator WOPPER for 
W pair production and decays into four fermions at high energy e + e~ colliders. 
The distinguishing feature of WOPPER is the inclusion of higher order electromag- 
netic corrections including soft photon exponentiation and explicit generation 
of exclusive hard photons. In contrast to fixed order calculations which have to 
be exponentiated by hand, WOPPER handles the multiphoton effects explicitly. 

The present version does not contain weak corrections. Forthcoming versions 
of the Monte Carlo generator will include weak corrections in the framework of 
effective Born cross sections Jss) . 
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A Distribution 

The latest release of WOPPER is available by anonymous ftp from 

crunch. ikp.physik.th-darmstadt .de 
in the directory 

pub/ ohl/ wopper 



23 



or on the World Wide Web at the URL 

\protect\vrule widthOpt\protect\href {http : // crunch. ikp .physik. th-darmstadt .de/monte-carlos . 

Important announcements (new versions, fatal bugs, etc.) will be made through 
the mailing list 

wopper-announce@crunch . ikp . physik . th-darmstadt . de 

Subscriptions can be obtained from 

majordomo@crunch . ikp .physik. th-darmstadt .de 

(send a message consisting of help to majordomo for instructions on how to 
subscribe, don't send such messages to the list itself). 

B Installation 

WOPPER is distributed in PATCHY format ||]. Plain FDRTRAN-77 versions can 
be made available on request. The source conforms to ANSI X3. 9-1978, with 
the exception of use of the following rather common extensions: 

• implicit none (may be disabled), 

• double complex arithmetic (essential), and 

• do ... end do constructs. 

Therefore, it should run without modifications on all platforms. 

B.l UNIX Systems 

On UNIX systems, the configuration, compilation and installation can be per- 
formed automatically according to the following sequence: 

$ ./configure 
$ make 

$ make install 

Optional interface code for external hadronization Monte Carlos, analysis rou- 
tines, etc. can be selected by giving command line options to the configure 
script. Figure |^ shows some of these command line options. This configure 
script has been created by the popular GNU Autoconf package and should 
work on all UNIX variants. 

If the — with-hepawk option has been specified in the configure step, a 
self test can be performed by make test. 

B.2 Non-UNIX Systems 

For non-UNIX systems configuration and compilation has to be performed man- 
ually from the CARDS file and the cradles. 
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Usage: configure [options] [host] 

Options: [defaults in brackets after descriptions] 
Configuration: 



— cache-f ile=FILE 

— help 

— no-create 
Host type : 

— build=BUILD 

— host=H0ST 

— target=TARGET 
Features and packages : 

— disable-FEATURE 

— enable-FEATURE [=ARG] 

— with-PACKAGE [=ARG] 

— without -PACKAGE 
— enable and 

— with-g77 



cache test results in FILE 

print this message 

do not create output files 

configure for building on BUILD [BUILD=HOST] 
configure for HOST [guessed] 
configure for TARGET [TARGET=HOST] 



do not include FEATURE (same as 
include FEATURE [ARG=yes] 
use PACKAGE [ARG=yes] 
do not use PACKAGE (same as 
with options recognized: 

use GNU Fortran 77 



-enable-FEATURE 



-with-PACKAGE=no) 



-enable-verbose-patchy display all patchy output 
-enable-internal 
-enable-not ime 
-enable-pedantic 
-with- 1 ibpath=P ATH 
-with-srcpath=PATH 
-with-hepawk 
-with-jetset 
-with-herwig 
-with-ariadne 
-with-circe 
-with-ranlux 
-with-cernlib 
-with-lepbench 
-enable-debug 
-enable-paper-a4 
-enable-paper-us 



do not use this! 
do not use timing functions 
no IMPLICIT NONE 
use PATH for libraries 
use PATH for source files (CARs) 
use HEPAWK for event analysis 
use JETSET hadronization 
use HERWIG hadronization 
use ARIADNE hadronization 
use CIRCE for beamstrahlung 
use RANLUX generator 
use CERNLIB 

link in the LEP2 benchmarking code 
compile for debugging 
use European (A4) paper 
use US (letter) paper 



Figure 8: Some of the comandline options of the configure script for WOPPER 
on UNIX systems. 

C External Symbols: Common Blocks and Sub- 
routines 

To avoid possible name clashes with other packages, all external symbols ex- 
ported by WOPPER proper begin with the two letters WW, except for the routine 
wopper itself and the /hepevt/ common block. 

• Common Blocks: 

The following common blocks are used by WOPPER: 

— /hepevt/, /hepspn/: standard common blocks for passing gener- 
ated events [Q 

— /wwpcom/: main parameter common block, holds all physical parame- 
ters. Application programs should access this common block through 
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the wwpsrv routine 

— /wwcbrn/ : internal parameters used for born cross section 

— /wwcdec/: internal parameters used for W decays 

— /wwcevt/: internal parameters used for event generation 

— /wwcsta/: statistics 

— /wwctri/: storage for keyword lookup 

Driver Program: 

— wwdriv: sample main program, which reads commands from stan- 
dard input and feeds them into wwdcmd 

— wwdloo: command loop, reading command from a terminal or file 
and executing them 

— wopper: the low level entry point into the generator for application 
programs 

— wwdcmd: WDPPER's command interpreter, the preferred entry point for 
application programs. Executes a single command 

— wwdlxi , wwdlxd, wwdlxs: Utility routines: tokcnization of input 

— wwdsig: UNIX signal handler 

— isatty: check if this job is run interactively 

Parameter Management: 

These routines are used to control the parameters common block /wwpcom/: 

— wwpsrv: server handling parameter changing requests 

— wwpini: block data supplying default values 

— wwpprn: print parameters 

Initialization: 

— wwinit: main entry point for initializations 

— wwigsw: initialization of electroweak parameters 

— wwicut: initialization of internal Monte Carlo parameters 

— wwibmx: find maximum of Born cross section 

— wwibmy: wrapper for wwuamo 

— wwibn: auxiliary function for finding maximum of on-shell cross sec- 
tion 

— wwibns: auxiliary function for finding maximum of cross section for 
symmetric virtualities 

— wwibna: auxiliary function for finding maximum of cross section for 
one W on-shell 

— wwibnv: auxiliary function for finding maximum of cross section for 
general case 

Final calculations: 

— wwclos: calculate total cross section, errors and close the generator 

— wwstat: statistics 

Hard Subprocess Generation: 
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— wwgen: main entry point for hard subprocess generation 

— wwgcfs: selection of initial and final state quantum numbers 

— wwgww: generation of angular distribution of pseudo W's 



— wwgppr 

— wwgdec 

— wwdqfl 



generate four-momenta of pseudo W's 

generate final state fermions from pseudo W decays 

select quark flavors in W decay 



Branching: 

— wwbini: generates the initial state photon radiation 
Cross Sections: 

— wwheli: coefficients of the helicity amplitude decomposition 

— wwxtot: total off-shell cross section 

— wwxint: auxiliary function for integration over pseudo W masses 

— wwxdif : differential W cross section 

— wwxdmx: maximum estimate of differential W cross section 

— wwxhel: helicity amplitude for e + e~ — ► 4/ 

— wwxhmx: maximum estimate of helicity amplitude for e + e~ — ► 4/ 
Accessing /hepevt/: 

— wweeni: enter identification of Monte Carlo and run 

— wweens: write summary record to /hepevt/ 

— wweent: enter one particle into /hepevt/ 

— wwenul: enter null particle into /hepevt/ 

— wwenew: clear /hepevt/ 

— wwecpy: copy a /hepevt/ entry 

— wwesf t: shift /hepevt/ entries 

— wwepsh: push /hepevt/ onto a stack (one level deep) 

— wwepop: pop /hepevt/ from a stack 

Hadronization: 

— wwpart: leave partons alone 

— wwlund: JETSET interface code 

— Iu4f rm: JETSET interface code 

— wwhwig: HERWIG interface code 

— hw4fdo, hw4f in, hw4fgo, hw4f cc, hw4f cs: HERWIG interface code 

— hwaend: HERWIG abnormal end routine 

— wwaria: ARIADNE interface code 

Random Numbers: 
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— wwrgen: returns a double precision uniform deviate 

— wwurng: (subroutine version) 

— wwrmz: random number generator RANMAR 

— wwrtst: test the portability of the random number generator 

— wwrtmz: test RANMAR 

— wwrkng: Knuth's random number generator 

— wwrknu: Knuth's random number generator (buffered version) 

— wwrkns: seed Knuth's random number generator 

— wwrknl: set luxury for Knuth's random number generator 

— wwrkni: initialize Knuth's random number generator 

— wwrknt: test Knuth's random number generator 

Utilities: 

— wwumsg: messages and error exit 

— wwulwr: convert input to lower case 

— wwuboo: boost a four vector 

— wwutim: still available CPU time for this job 

— wwuamo: multidimensional minimization 

— wwumin: onedimensional minimization 

— wwuons: Gram-Schmidt procedure 

— wwuort: another Gram-Schmidt procedure 

Canonical cuts: 

— adloth: apply canonical cuts in /hepevt/ 

— adloip: inner product of vectors in /hepevt/ 

— adloan: angle to beam in /hepevt/ 

— adload: add four momenta in /hepevt/ 

LEP2 benchmarks ||: 

(optional, selected with configure — with-lepbench) 

— hepawk: main analysis routine 

— lbf ill: fill a histogram 

— lbcheb: Checbycheff moments 

— lbpowr: power moments 

— lbprin: print output 

— lblsct: select event 

Keyword search: 
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(using the dynamic tries described in [ETJ) 

— wwtins: insert a new keyword 

— wwtlup: look up a (possibly abbreviated) keyword 

— wwtnew: insert new a node into the trie 

— wwtlen: calculate length of keyword 

— wwtc2a: convert keyword from character* (*) to integer (*) 
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Test Run 

WOPPER version 1.5 is distributed together with a sample command file and 
hepawk script, which are given below. To run this example, the user will need to 
link WOPPER with the CERN library, because histogramming is done by HBOOK 



The file sample . wopper is read from standard input (unit stdin, which is 
initialized to 5) , and sample . hepawk is read from the file SCRIPT (i.e. under MVS 
from the file which has been allocated to the DDNAME SCRIPT and under UNIX 
from the file script or from the value of the environment variable SCRIPT). 

sample .wopper 

Here is a simple WOPPER command file, setting up parameters and generating 
10000 events. 

# sample .wopper — sample WOPPER command file 

# parameters 
set ebeam 87.5 

# run 
init 

gen 10000 

close 

quit 

sample .hepawk 

This is a small hepawk |2?| script that counts the muons from the W decays 
and plots a histogram of their energy distribution. The first generated event is 
dumped to illustrate the usage of the /hepevt/ common block. 

# sample .hepawk — sample HEPAWK analyzer for WOPPER. 



[ID- 



BEGIN { 



printf ( 

printf ( 

printf ( 

printf ( 

printf ( 



"\nWelcome to the WOPPER test:\n"); 
"Monte Carlo Version: °/.s\n" , REV); 



Run: y„d\n", RUN); 
Date : 7„s\n\n" , DATE) ; 
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E_max = 100; 
N_chan = 50; 
h_muon_ energy 

= bookl (0, "Muon-energy" , N_chan, 0, E_beam) ; 
nr_muons =0; # initialize counter 

dumped_an_event = 0; 



{ 

if (dumped_an_event == 0) { 

dump ("vs") ; # Dumping first event 

dumped_an_event++ ; 

} 



for (@p in LEPTQNS) 

if (abs(@p:id) == _pdg_muon) { 
fill (h_muon_energy , @p:p:E); 
nr_muons++ ; 

} 

} 



END { # Dump some numbers 

printf ("\nRESULTS:\n********\n\n") ; 

printf ("Total events: °/,d, total cross section: °/,g pb\n" , 

NEVENT, XSECT * le9) ; 
printf ("Number of muons : °/,d\n\n" , nr_muons) ; 
plotO ; 

printf ("\ndone.\n") ; 
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sample . output 

The following output should result from the input files above, up to small round- 
off errors and different FDRTRAN-77 default output formats. 



message 
message 
message 
message 
message 
message 
message 
message 
message 
message 
message 
message 
message 
message 
message 
message 
message 
message 
message 
message 
message 
message 
message 
message 
message 
message 
message 
message 
message 
message 
message 
message 
message 
message 
message 
message 
message 
message 



Starting WOPPER, Version 1.05/00, (build 960530/2055) 
linked with ARIADNE, 
linked with CIRCE, 
linked with HEPAWK. 
linked with HERWIG. 
linked with JETSET. 
enabled RANLUX. 
starting HEPAWK, Version 1.6 

"GF scheme" selected: 

Using GFERMI and ALPHA as input, calculating SIN2W. 

Using derived W width. 

Parameters used in this run: 

AHPLA = 128.00000 ( = l/alpha(2 M_W) ) 

SIN2W = 0.23098 (effective mixing angle) 

GFERMI = 0.116639E-04 GeV**(-2) 

ALPHAS = 0.12000 

GAMM1W = 2.08780 GeV (S.M. value, used) 

CKMVUS = 0.21960 CKM Matrix (Cabibbo angle) 

CKMVCB = 0.04000 CKM Matrix 

CKMVUB = 0.00320 CKM Matrix 

Z-e-e couplings : 

g_V = -0.01414 
g_A = -0.18586 
W-e-nu coupling: 

g = 0.23050 
Z-W-W coupling: 

g_ZWW = 0.57171 
gamma-W-W coupling: 
g_gWW = 0.31333 

Using standard model triple gauge couplings. 
Using energy-dependent width for W and Z propagators. 
Including "CCll"-type background diagrams. 
Including initial state radiation (leading logs) . 
No Coulomb correction. 

Selected parton level events. 



Welcome to the WOPPER test: 



Monte Carlo Version: v01.05 (May 30 00:00:00 
Run: 1035996352 
Date: May 30 20:57:00 1996 



1996) 



Dumping the event record for event # 
There are 13 entries in this record: 
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Entry # 1 is an incoming (HERWIG convention) electron 

p: (0.8750E+02; 0.0000E+00, 0.0000E+00, 0.8750E+02), m: 0.0000E+00 



Entry # 2 is an incoming (HERWIG convention) positron 
p: (0.8750E+02; 0.0000E+00, 0.0000E+00, -.8750E+02), m: 





. 0000E+00 


Entry # 3 is the CMS system (HERWIG convention) 
p: (0.1750E+03; 0.0000E+00, 0.0000E+00, 0.0000E+00), 


m: 





. 1750E+03 


Entry # 


4 


is a null entry. 








Entry # 


5 


is a null entry. 








Entry # 


6 


is a null entry. 








Entry # 


7 


is reserved for model builders. 








Entry # 


8 


is reserved for model builders. 








Entry # 9 is an existing photon 

p: (0.1554E+00; 0.6972E-03, - . 1679E-03 , -.1554E+00), 
The mother is the positron # 5. 


m: 


0. 


. 0000E+00 


Entry # 


10 


is an existing muon-neutrino 









p: (0.5820E+02; 0.5040E+02, 0.2417E+02, -.1620E+02), m: 0.0000E+00 
The first mother is the W+ #7. 
The other mother is the anti-muon # 11. 



Entry # 11 is an existing anti-muon 

p: (0.2788E+02; -.2410E+02, -.6728E+01, 0.1229E+02), m: 0.1057E+00 
The first mother is the W+ #7. 
The other mother is the muon-neutrino # 10. 



Entry # 12 is an existing muon 

p: (0.4687E+02; -.2813E+00, -.3362E+02, 0.3265E+02), m: 0.1057E+00 
The first mother is the W- #8. 
The other mother is the anti-muon-neutrino # 13. 



Entry # 13 is an existing anti-muon-neutrino 

p: (0.4191E+02; -.2602E+02, 0.1617E+02, -.2859E+02), m: 0.0000E+00 
The first mother is the W- #8. 
The other mother is the muon # 12. 



RESULTS : 

Total events: 10000, total cross section: 13.06 pb 

Number of muons : 2195 

lMuon-energy 

HB0QK ID = 1 DATE 30/05/96 
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172 2 

168 X 

164 X 

160 7 X 

156 X X 

152 X5 X5 

148 XX XX 

144 XX2XX 

140 2 XXXXX 

136 7X7XXXXX 

132 XXXXXXXX 

128 XXXXXXXX 

124 XXXXXXXX7 

120 XXXXXXXXX 

116 XXXXXXXXX 

112 2 XXXXXXXXXX 

108 X 7XXXXXXXXXX 

104 5X XXXXXXXXXXX 

100 XX XXXXXXXXXXX 

96 7 XXXXXXXXXXXXXX 

92 X XXXXXXXXXXXXXX 

88 X XXXXXXXXXXXXXX 

84 X XXXXXXXXXXXXXXX 

80 X XXXXXXXXXXXXXXX 

76 X2XXXXXXXXXXXXXXX 

72 XXXXXXXXXXXXXXXXX 

68 XXXXXXXXXXXXXXXXX 

64 XXXXXXXXXXXXXXXXX 

60 XXXXXXXXXXXXXXXXX 

56 XXXXXXXXXXXXXXXXX 

52 XXXXXXXXXXXXXXXXX 

48 XXXXXXXXXXXXXXXXX 

44 XXXXXXXXXXXXXXXXX 

40 XXXXXXXXXXXXXXXXX 

36 XXXXXXXXXXXXXXXXX 

32 XXXXXXXXXXXXXXXXXX5 

28 XXXXXXXXXXXXXXXXXXX 

24 XXXXXXXXXXXXXXXXXXX 

20 2XXXXXXXXXXXXXXXXXXX 

16 xxxxxxxxxxxxxxxxxxxx 

12 XXXXXXXXXXXXXXXXXXXX72 

8 XXXXXXXXXXXXXXXXXXXXXX 
4 25 7XXXXXXXXXXXXXXXXXXXXXXXX7 5 

CHANNELS 10 1 2 3 4 5 

1 12345678901234567890123456789012345678901234567890 

CONTENTS 100 11 11111111111 

10 139700901333554652831 
1. 12 34725329672575901903401943 2 

LOW-EDGE 10 111112222233333444445555566666777778888899999 
1 . 2468024680246802468024680246802468024680246802468 
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* ENTRIES = 2195 * ALL CHANNELS = 0.2195E+04 * UNDERFLOW = 0.00 

* BIN WID = 0.2000E+01 * MEAN VALUE = 0.4615E+02 * R . M . S = 0.10 



done. 

wwdriv: message: bye. 

D Revision History 

Version 1.5, June 1996 

• Anomalous couplings. 

• "CC11" diagrams. 

• LEP2 standardized ARIADNE, JETSET and HERWIG interfaces. 

Version 1.4, Fall 1995 

• LEP2 workshop. 

Version 1.3, April 1995 

• Canonical cuts and input parameters. 

• Fixed inconsistent phase conventions, which resulted in wrong angular 
distributions. 

Version 1.2, July 1994 

• Coulomb correction. 

• Improved JETSET and HERWIG interfaces. 

Version 1.1, February 1994 

• Hadronization, JETSET and HERWIG interfaces. 

• Minor bug fixes. 

Version 1.0, 1993 

First official release, submitted to the Computer Physics Communication Li- 
brary. 
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